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Abstract 

In this paper we investigate the vacuum polarization effects associated with quantum 
fcrmionic charged fields in a generalized (d+ 1)— dimensional cosmic string space-times con- 
sidering the presence of a magnetic flux along the string. In order to develop this analysis 
we calculate a general expression for the respective Green function, valid for several different 
values of d, which is expressed in terms of a bispinor associated with the square of the Dirac 
, operator. Adopting this result, we explicitly calculate the renormalized vacuum expectation 

' values of the energy- momentum tensors, (T^)fl en ., associated with massless fields. Moreover, 

for specific values of the parameters which codify the cosmic string and the fractional part of 
the ratio of the magnetic flux by the quantum one, we were able to present in closed forms 
the bispinor and the respective Green function for massive fields. 
PACS numbers: 98.80. Cq, ll.10.G7i, 11.27. + d 
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1 Introduction 



^ . Braneworld model has received renewed interest in recent years. By this scenario our world is 

represented by a four dimensional sub-manifold, a three-brane, embedded in a higher dimen- 
sional spacetime [H [2] . The idea that our Universe may have more than four dimensions was 
proposed by Kaluza [3J, with the objective to unify gauge theories with gravitation in a geomet- 
ric formalism. Of particular interest are the models introduced by Randall and Sundrum [H [5] . 
The corresponding spacetime contains two (RSI), respectively one (RSII), Ricci-flat brane(s) 
embedded on a five-dimensional Anti-de Sitter (AdS) bulk. It is assumed that all matter fields 
are confined on the branes and gravity only propagates in the five dimensional bulk. In the RSI 
model, the hierarchy problem between the Planck scale and the electroweak one is solved if the 
distance between the two branes is about 37 times the AdS radius. 

Although topological defects have been first analysed in four-dimensional spacetime [3 [7], 
they have been considered in the context of braneworld. In this scenario the defects live in 
a n— dimensions submanifold embedded in a (4 + re)— dimensional Universe. The domain wall 
case, with a single extra dimension, has been considered in [2j. More recently the cosmic string 
case, with two additional extra dimensions, has been analysed [HI [9] . It has been shown that 
the gravitational effects of global strings can be responsible for compactification from six to 
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four spacetime dimensions, naturally producing the observed hierarchy between electroweak 
and gravitational forces. 

Cosmic strings are topologically stable gravitational defects which appear in the framework 
of grand unified theories. These objects could be produced in very early Universe as a result 
of spontaneous breakdown of gauge symmetry [6] [TJ. Although the recent observation data 
on the cosmic microwave background have ruled out cosmic strings as the primary source for 
primordial density perturbation, they are still candidate for the generation of a number of 
interesting physical effects such as gamma ray burst [lOj . gravitational waves [11] and high 
energy cosmic ray [12]. Recently, cosmic strings have attracted renewed interest partly because 
a variant of their formation mechanism is proposed in the framework of brane inflation |13j-|15j. 

The simplest theoretical model describing an idealized cosmic string, i.e., straight and in- 
finitely thin, is given by a delta-type distribution for the energy-momentum tensor along the 
linear defect. As the solution of the Einstein equation, the geometry of the space-time produced 
by this source presents a conical singularity for the curvature tensor on its top. Under classical 
field theory viewpoint, this object can also be formed coupling the energy-momentum tensor 
associated with the Higgs U(l)— gauge system investigated by Nielsen and Olesen [16J with the 
Einstein equations. This project was successfully analysed by Garfinkle [17] . He found static 
cylindrically symmetric solutions representing vortices, as in flat space-time, and shown that 
asymptotically the space-time around the vortices is a Minkowski one minus a wedge. Their 
core have a non-zero thickness, and the magnetic fields vanishes outside them. Two years later 
Linet [18] obtained, as a limit case, exact solutions for the metric tensor and Higgs field. He was 
able to show that the structure of the respective space-time corresponds to a conical one, with 
the conicity parameter being expressed in terms of the energy per unity length of the vortex. 

The vacuum polarization effects associated with scalar and fermionic fields in the space-time 
of an idealized cosmic string, have been analyzed in [19]- [24], and [25]-[28], respectively. It has 
been shown that these effects depend on the parameter which codify the conical structure of the 
geometry. Moreover, considering the presence of magnetic flux along the cosmic strings, there 
appears an additional contributions to the vacuum polarization effect associated with charged 
fields [29]- [33]. Another type of vacuum polarization takes place when boundaries are presents. 
The imposed boundary conditions on quantum fields modifies the zero-point energy fluctuations 
spectrum and result in additional shifts in the vacuum expectation values of physical quantities, 
such as energy density and stress. This is the well-known Casimir effectslJ The analysis of 
Casimir effects in the idealized cosmic string space-time have been developed for a scalar [35] 
and vector fields [36], obeying boundary conditions on the cylindrical surfaces!! 

The investigation of quantum effects in corresponding braneworld models is of considerable 
phenomenological interest, both in particle physics and cosmology. Quantum effects provide 
natural alternative for stabilizing the radion fields in a braneworld. The corresponding vacuum 
energy gives contribution to both the brane and bulk cosmological constant and, hence, has 
to be taking into account in the self-consistent formulation of the braneworld dynamics. Re- 
cently the fluxes by gauge fields play an important role in higher dimensional models including 
braneworld scenario (see for example [38] ) . They provide stabilization mechanism for all moduli 
fields appearing, in particular, in various string compactifications. Motivated by these facts, we 
decided to analyse the vacuum polarization effects associated with a charged spin— 1/2 field in a 
higher-dimensional cosmic string space-time considering an infinitely thin magnetic flux running 
along it. Specifically we are interested to calculate the vacuum expectation value (VEV) of the 
energy-momentum tensor. Although, as we shall see, for a {d + 1)— dimensional cosmic string 
space-time, the renormalized VEV of the energy-momentum tensor associated with massless 

For a review, see Ref. [34] 
2 Also vacuum polarization effects induced by a composite topological defect has been analysed in [37] 



2 



quantum fields presents the general form below 





(1) 



r d+l ' 



being r the azimuthal distance to the string, the components of the tensor Fg have different 
expressions depending on the dimension of the manifold. 

This paper is organized as follows: In Section[2]we present a general expression for a charged 
massless spin— 1/2 Green function valid for an arbitrary dimensional cosmic string space- 

time, d>2. Although our main objective is to investigate the vacuum polarization effects in a 
six-dimensional spaces, the Green function obtained allows us to develop this analysis for lower 
dimensions as well. Moreover, in this section we also present a general expression for Green 
function associated with massive field for a particular choice of the parameters which codify 
the presence of cosmic string, and the fractional part of the ratio of the magnetic flux by the 
quantum one. In this way, the massive fermionic Green function is expressed in terms of a finite 
sum of modified Bessel function K v , allowing us to present the renormalized VEV of the energy- 
momentum tensor in a closed form. In the Section [3] we investigate the spinor Green functions 
in coincidence limit and extract all divergences in manifest form. Using this result, we provide 
explicit expressions for all components of the renormalized VEV of the energy-momentum tensor 
for different dimensions of the space, and for massless and massive fields. In Section HJ we 
summarize the most important results obtained. In Appendix [A] we present a generalization 
of the generation function for the modified Bessel functions, needed to construct the fermionic 
Green function as a finite sum of images of the Minkowski space-time functions. The Appendix 
IB1 contains some technical formulas related with the calculation of the renormalized VEV of the 
energy-momentum tensor in Section [3l 

2 The spinor Green function 

This section is mainly devoted to calculate the Feynman propagator associated with a charged 
fermionic field propagating in a (d + 1)— dimensional cosmic string space-time whose geometry 
is given by the line element 



Specifically for a six-dimensional space, the coordinate system reads: x = (t, r, <p, x, y, z), with 
r > 0, (p G [0, 2tt], and t, x % G (— oo, oo), i = 3, 4, 5. a is parameter smaller than unity 
associated with the linear mass density of the string. In the braneworld scenario the space-time 
given by ([2]), the bulk, represents a conical 3— brane transverse to a two flat space. 

In order to develop the calculation of the spinor Green function we shall adopt the 8x8 
Dirac matrices T M given below, which can be constructed in terms of the 4x4 ones [39], [40|: 



where 75 = Z7°7 1 7 2 7 3 , I represents the 4x4 identity matrix and r and (p stand the ordinary unit 
vectors in cylindrical coordinates. This set of matrices satisfies the Clifford algebra {T M , F N } = 




(2) 






(8)- 
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In the analysis of the vacuum polarization effects, we also consider the presence of an ex- 
tra magnetic field running along the string. This magnetic field configuration is given by the 
following six-vector potential 



A M = Ad M <f (4) 



being A — ! ' 



2tt- 



The spinor Feynman propagator, defined as [H] 

iS F (x,x') = (0\T(^(x)^(x'))\0) , (5) 
with ^ = \jdr°, satisfies the non-homogeneous linear differential equation, 

{if + e4- M) S F (x, x') = 7^=^ " ^0^(8) , (6) 

where g = det(gMN)- The covariant derivative operator reads 

W = T M (d M + U M ) , (7) 
being Um the the spin connection given in terms of the T— matrices by 

n M = ~r N v M r N , (8) 

and 

4 = r M A M . (9) 

The Green function given in ([6]) is a bispinor, i.e., it transforms as ^ at x and as at x' . 
In [33] we have shown that if a bispinor, T> F (x,x'), satisfies the differential equation 



[□ - ieg MN (D M A N ) + ieZ MN F MN - 2ieg MN A M V N 
-e 2 g MN A M A N -M 2 -\n 



V F {x,x') = —^6 6 (x-x')I {8) , (10) 



with 

1Z being the scalar curvature and the generalized d'Alembertian operator given by 

□ = 9 MN V M V N = g MN {d M V N + n M Viv - {f/iv}Vs) , 

then the spinor Feynman propagator may be written as 

S F (x, x) = (if + + M) V F (x, x) . (12) 

Applying this formalism for the system under investigation the operator /C, which acts on 
the left hand side of pUj) . reads 



K = A+ ^( 1 -«) S (8)^-i^( 1 -«) 2 + ^( 1 -«)^(8) 

2ie . „ e 



- —j-^Adp - -x-zA* - M z , (13) 
where 

E ?8) = ( S 3 E 3 ),wi.hE3=(; 3 °) (14) 
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anc|l 

A = -d 2 + d 2 + ^d r + ^d 2 +d 2 x + d 2 + d 2 z . (15) 

Due to the fact of EjL be a diagonal matrix, the bispinor T>p(x,x') is diagonal, too. In this 
way we can obtain an expression for this Green function analyzing only the effective 2x2 matrix 
differential equation below: 



2ie Ad u> --^A 2 -M 2 



a 2 r 2 f o^r 2 



vP(x,x') = —^=6 6 (x-x')I (2) . (16) 



-g 

So the complete Green function is given in terms of T> F (x, x), in a diagonal matrix form. 

On basis of these results, for this six-dimensional cosmic string space-time, the fermionic 
propagator reads: 



Sp(x, x') 



iT°d t + iT r d r + iT*d v + iT% - l^T^f 8) + ^T* + M 



D F (x,x') . (17) 



The vacuum average value for the energy-momentum tensor can be expressed in terms of 
the Euclidean Green function. It is related with the ordinary Feynman Green function [41] by 
the relation T>E{r,r;r' ,r') = —iDf(x,x'), where t = it. In the following we shall consider the 
Euclidean Green function. 

In order to obtain the Euclidean Green function T> E (x,x') in explicit form, let us find the 
complete set of bispinor which obey the eigenvalue equation 

iC^ x {x) = -A 2 $ A (x) (18) 

with A 2 > and being K, the Euclidean version of the differential operator given in (|16p . So we 
may write 



V?M = E = f *EMx)»jM e-* ■ (1») 



The eigenfunctions of (fTKl) can be specified by a set of quantum number associated with operators 
that commute with K, and among themselves: p T = —id T , pi = —idi for i = 3, 4, 5 , L v = —id v 
and the spin operator, 173. Let us denote these quantum numbers by (k T , k l , n, a), where 
(k T , k l ) £ (—00, 00), n = 0, ±1, ±2, ... , a = ±1. Moreover, these functions also depend on 
the number p, which satisfies the relation A 2 = p 2 + k 2 + M 2 and assumes values in the interval 
[0, oo). 

Although we have developed this formalism for a six-dimensional cosmic string space-time, 
it can be adapted for a three, four and five dimensional spaces. The reason resides in the 
representations adopted for the Dirac matrices in these dimensions. For three dimensional case 
the Dirac matrices are given in terms of 2 x 2 Pauli matrices. For four and five dimensions, the 
Dirac matrices are given by the 4x4 off-diagonal matrices from the first four and five matrices 
given in ([3]). Of course for these cases we have to discard the derivative associated with the 
respective extra coordinates in (|15p . So on basis on these arguments it is possible to generalize 
the eigenfunction of the operator (fT8~j) by: 



s For this geometry the only non- vanishing spin connection is 11^ = i(l — a)E 



(8)- 
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^ = ^(^^i/a ^Vl/o^)^ (20) 



where 



" <+ K ° ) ' "'M ? ) < 21 > 

are the eigenfunctions of the the operator 03, J^(z) is the Bessel function, and v^- = n± — — 
[N + 7), where we have defined eA = e-^ = N + 7, the ratio of the magnetic flux by the 
quantum one, in terms of an integer number, N, and its fractional part, 7. In (|20p we may 
assume N = 1, 2, 3, 4, which are related with the dimensions of the space considered: three, 
four, five or six, as explained before. 

Now we are in position to calculate the Green function by using (fl~9l) , (f20l) and (j2Tj) as show 
below 

di ag (J ll/+l/a (pr)J ll/+l/a (pr'), J W -\ /a (pr)J w -\ /a (jrr')) e -( P 2 +fc 2 +*' 2 ) . ( 22 ) 
With the help of [32] we can express ([22]) by 



D E^) ~ a{4n) N/2 + l J s N/2 + i e Z. e 

diag(J K | /a (rr72a), I ]v -\ /a (rr' /2s)) , (23) 

where Iu(z) is the modified Bessel function. 

The calculations of the renormalized VEV of the energy-momentum tensor associated with 
massive quantum fields in cosmic string space-time, have been developed by many authors cited 
in |25j, 129] . All of these calculations involve very complicated integrals. In order to present 
more suitable expressions, the authors make applications taking the massless limit in previous 
results. Here, in this paper, we shall explicitly exhibit the renormalized vacuum expectation 
value of the fermionic energy-momentum tensor for different dimensions of the space-time in 
closed expressions. So, in order to do that we shall consider two different limiting cases which 
allow us to reach this objective as shown in the following subsections. 

2.1 Massless fermionic fields 

In this subsection we shall calculate the bispinor associated with a charged massless fermionic 
fields in an arbitrary dimensional cosmic string space-time. So taking M = in (|23|) we obtain 



V%\x,x') 



N+3 N N—l 

a(2ir)^~ (rr')~ (— sinhw)^ - 

JV-l 



^ e ^-v)[ %h«-i/2^ u ) ^ 1 (24) 

Q\u-\/a-l/2( COshu )) 
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where 

C o S ^u= iAx)2 + r * + r ' 2 . (25) 
2rr' y ' 

In the above expression, Q x is the associated Legendre function. It is possible to express this 
function in terms of hypergeometric function [12] by: 

r( u _i_ \ _i_ i p -(v+\+i)u 
Q x (coshu) = e iXn 2 x y^F y\ + - _ 9 - — j- (sinh u) x x 



r(^ + 3/2) (i_ e -2«)A+V2 

1 



F(A + l/2, -A + l/2; I/ + 3/2; ^ ) . (26) 



In this analysis we identify the parameter f with Iz/^l/a — 1/2 and take A = (N — l)/2. So the 
relevant hypergeometric function is 

„(N N - 2 . , . . 1 
F I — , — ; |^|/a + l ; 



2 ' 2 ' 1 " ' 1 - e 2u . 

If iV is a even number, this function becomes a polynomial of degree however being N an 
odd number, this function is an infinite series. 

As special applications of this formalism, let us consider specific values for N as shown in 
the following. 

2.1.1 N=l 

For this case the manifold corresponds to a three-dimensional cosmic string space-time. The 
Green function can be obtained by substituting the Legendre function by its integral represen- 
tation below 1 42 1, 



1 f°° e~ 
Qu- 1/2 (cosh u) = —= / dt , 
ylJu vcosni — cosnu 



(27) 



in (|24p and developing the sum over n. After some steps we get 

(2) e iN(<P~<p') ( f°°dt . T )(t) C 
JY ■ 'It' t\ — —— I Vcosht-coshu 

v ' \ J « Vcosht-coshu 



(28) 



where 



S (±) (t) = gg^j sinh(^t/q) - sinh[(^ - l)t/a] ^ 
cosh(i/a) — cos(ip — <//) 

being 6^ = » =F 7- 
2.1.2 N=2 

For the case with N = 2 the bulk corresponds to a four-dimensional cosmic string space-time. 
For this case the Green function is very well known. It takes a simple form due to the fact that 
the respective associated Legendre functions be expressed by 

Qy_ 2 1/2 (cosh.)=^-|^=. (30) 
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Substituting this function into ([21]) it is possible to develop the sum on the angular quantum 
number n, providing 

.(2), , e <*W) /£(+)(„) \ 

£lX ' XJ S^Wsinh^V £<->(«)) / ' 

(31) 

with S^-^w) having the same functional form as in (1291) , 

2.1.3 N=3 

This case is a new one, it corresponds to a five dimensional cosmic string space-time having a 
magnetic flux running along the string. For this case the bispinor depends on the associated 
Legendre function Q^(z), which can be expressed in a integral representation by using the 
relation Ql(z) = {z 2 — i) 1 / 2 and ([27|) . After some intermediates steps we obtain 



s<--Ht) 

Vcosh t—z 



=cosh u 



(32) 



2.1.4 N=4 



This case is also a new one, it corresponds to a six-dimensional cosmic string space-time. The 
complete Green function is expressed by a 8 x 8 diagonal matrix. As in (|2.1.2p , this function can 

(2) 

also be written in terms of elementary functions. The reason is because for N = 2, T> E '{x' ,x) 
depends on Q^i^C 2 )' which is given by 



„3/2 , , s . / vr (cosh « + sinh u) . . 

Q_ 1/2 (coshn) = -^^- i -e _ . (33) 

Substituting this expression into (|24p . and after some intermediate steps we arrive to: 



Vf(x',x) = * N0 ^x 

1671^ a(rr') 2 sinh u 



cosh uS( + \u) — sinh -uS 1 '^ (u) 

cosh uS^ (u) — sinh S'^ (u) 

(34) 

where the prime denotes derivative with respect to u. 

2.2 Special case where a = - being q an integer number 

In this subsection we shall provide a closed expression for the bispinor considering a non- 
vanishing mass term for the fermionic field. This will happen for a very special situation when 
a is equal to the inverse of an integer number. Smith, Davies and Sahni, and Souradeep and 
Sahni, cited in [19], have shown that when this occur, i.e., for a = |, the scalar Green functions 
can be given as a sum of q images of the Minkowsiki space-time functions. Recently the image 
method was also used in [35] to provide closed expressions for massive scalar Green functions for 
a higher-dimensional cosmic string space-time. The mathematical reason for the use of image 
method in these applications resides in the fact that the order of the Bessel functions which 



S 



appear in the derivations of the Green functions, become an integer number. Unfortunately, 
for the fermionic case, the order of the Bessel function depends, besides the integer angular 
quantum number n, also on the the factor ^ 7 ■ However, if we consider a charged fermionic 
field in the presence of a magnetic flux running along the string, an additional effective extra 
factor will be present, the fractional part of the ratio of the magnetic flux by the quantum one, 
7. For the case where 7 is equal to — , and a = 1/q, the order of the Bessel function be- 
comes an integer number, and, in this case it is possible to use the image method to obtain the 
fermionic Green function in a closed form. Although being a very special situation, the analysis 
of vacuum polarization effects in this circumstance may shed light on the qualitative behavior 
of these quantities for non-integer q. 



2.2.1 Heat kernel 

Let us now investigate the effective 2x2 diagonal matrix heat kernel, given in the integrand of 
(US]). Using fl23|), we have 

Defining 

KM(x,x';s) = Y,e mAip I W+{/a (rr'/2 S ) (36) 

n 

with v + = n + — — N — 7, we can see that for 7 = — and assuming a = |, the order of 
the modified Bessel function becomes an integer number if q is also an integer number. In this 
way (f36|) can be expressed by [43jFl 

KW(x, x'; s) = e lNA * £ e inA *I nq (rr' /2s) = e iR ^ - £ cos (^+¥) . (37) 

n 1 k=0 

Analogously, defining 

K(-\x,x';s) = Y,e in ^I\ v -\ /oi {rr'/2s) (38) 

n 

with v~ = n — ^ — N — 7 = n — N — 1 + a, we obtain 

K(-\x,x';s) = e^+^^e^Inq+xirr'^s) . (39) 

n 

Using the recurrence relation, I u+ i(z) = I' u (z) — jl v {z) and Eq. (f3T|) . we get, after some steps, 
the following expression: 

K(-\x,x';s) = e Wi-V,)A, lg e ^-(^+¥) e -^ . (40 ) 

^ fc=0 

The Green function is given by 

poo 

T>%\x',x)= ds K(x,x';s) . (41) 



4 See the proof in Appendix El 
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Using (|57|) . ([1DD and with help of [32], this function becomes 



^V, *) = (o^)W M7V/2 E wm K wWPk) [ ei (i-i /g) A^- 2 ^/ g J , (42) 



iNA(p I' 1 
fc=0 

with 

p 2 . = (Ax) 2 + r 2 + r /2 - 2rr' cos (Ay>/g + 2-irk/q) , (43) 

being if^ the modified Bessel function. 

As we can see, the above Green function of the square of the Dirac operator corresponds, 
up to a phase factor, to a sum of q images of Minkowski space-time functions. For q = 1, and 
N = 2, the above expression reduces to the standard four-dimensional Green function 

V%\x>, x) = -^-- Kl (M Po ) J (2) . (44) 

We want to finish this section by saying that having obtained the complete Green functions 
associated with the square of the Dirac operator, for massless and massive fields, the fermionic 
propagators can be given by applying the Dirac operator on this bispinor according to (fl~2|) . 



3 Vacuum expectation value of the energy-momentum tensor 

In this section we shall calculate in a explicit form the renormalized vacuum expectation value 
of the energy-momentum tensor, {Tg) Ben- . As we have already mentioned in the last section, 
two physical situations will be considered here. Because the metric tensor does not depend of 
any dimensional parameter, for the massless fermionic field case we can infer that the VEV of 
the energy-momentum tensor will depend only on the radial distance r. For massive fermionic 
case, it is expected a dependence on the mass of the field too. 

The renormalized VEV of the energy-momentum tensor must obey the conservation condi- 
tion, 

y A {T^) Ren . = , (45) 
and for massless fields provide the correct trace anomaly for space-time of even dimension [44] : 

(Tjt)Ren. = ^rzrTr ai+d(x,x) . (46) 

(4tt) 2 2 

However, because this space-time is locally flat and there is no magnetic field in the region 
outside the string, the coefficients 02 and 03 vanish [35j [38] . El As we shall see, taking into 
account these informations, it is possible to express all components of the energy-momentum 
tensor in terms of only one. Here we shall calculate the energy density only. 

Using the point-splitting procedure [41] , the VEV of the energy-momentum tensor has the 
following form: 

(T AB (x)) = \ lim Tr [T A (D B - D B .) + F B (D A - D A ,)] S F (x,x') , (47) 

5 For massive fields the trace of the VEV of energy-momentum tensor does not vanish. In fact it is (T^)uen. = 
M 2 Tr T> Ren . (x, x). 
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where Dm = Vm — ieAw, and the bar denotes complex conjugate. Because the dependence 
of the fermionic Green function on the time variables, the zero-zero component of the energy- 
momentum tensor reads: 



(Too(x)} = lim Tr T d S F (x, x) 



(48) 



which can be expressed by 



(T 00 (x)) = -i lim d 2 Tr V F (x',x) 



lim d 2 Tr V E (x,x') . 

r 1 — 



(49) 



X 



X 



In the obtainment of the above expression we have to use the fact that the bispinor T>f{x' , x) is 
diagonal and Tr r°PX> E (:r, x') = 0. 

Now after these brief comments about general properties of the VEV of the energy-momentum 
tensor, we shall start explicit calculations of this quantity for massless and massive fermionic 
fields specifying the dimension of the manifold. 

3.1 Massless case 

In this subsection we shall calculate the zero-zero component of VEV o the energy-momentum 
tensor for massless fermionic fields, specifying the dimension of the manifold. Let us start with 
a three dimensional cosmic string space-time. 

3.1.1 N=l 

The vacuum polarization effects associated with massive scalar fields in a three-dimensional con- 
ical space-time has bee developed by Guimaraes and Linet cited in [19]. Because the calculation 
of the VEV of the energy-momentum tensor involves complicated integrals, the authors give its 
expression only for the massless case, which coincides with with the result found by Souradeep 
and Sahni, also cited in [19]. However, as far as we know, no calculations regarding to fermionic 
fields in the presence of a point-like magnetic field have developed in this space-time. 
For this three-dimensional cosmic string space-time the Dirac matrices reads: 



where cosh-u = +r ■ Let us take first r' = r, Aip = in S^\t) given by ([29]) . and define 

new function, 



7° = cr 3 , 7 1 = if ■ a and j 2 = — <2> ■ a , (50) 

ar 

being a k the Pauli matrices. In this space-time the Euclidean bispinor, Ve(x',x), is the 2x2 
diagonal matrix given below 




(51) 




(52) 




(53) 
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being z = cosh u = 1 + . 

In order to calculate the VEV of the zero-zero component of the energy-momentum tensor, 
let us write the function I^\z) as: 

JW(,) = l[ ± \ z)+ t\ z )=2[ 1 dy 



S<±>(y) 



v / 2TFv / y 2 -(^- 1 ) 



+ 2 



£(±)(y) 



V2T^Vy 2 -(^-i) ' 



(54) 



where S^(y) = S^ ± \t(y)). Subtracting and adding into the integrand of l[ the two first 
terms of the power expansion 



_a | [(1- a 2 ) +6^-1)] ^ | 



V2 + 2/ 2 y 



6a 



we get 



r(±) _ r(±)/*« , A±)sing 
1 1 — J l "T J l 



where 



7 (±)/- (z) =2 



riy 



and 



J 



(±)sing 
1 



(*) = 2 



t Vy 2 - (* - i) 

dy 



<g (±) (y) _« 
>/2 + y 2 y 



[(l-a 2 ) + 6«5 ± (<5 ± -l)] 



6a 



a [(l-a 2 ) + 6^(^ -1)] ' 
y 6a 



(55) 



(56) 



(57) 



(58) 



As we shall show, I^f in together with 1% \ provide finite contributions to the VEV. The 



only divergent contribution is contained in j^ mn9 ) which has the form 



I 



(±)sing 
1 



2a(3 



+ V2 



■ [(l-a 2 ) + 65 ± (5 ± -1)} 
3a 



(59) 



being 



(3 = arctg 



2-z 
z-1 



and z = 1 + 



(At) 
2r 2 



(60) 



Expanding jj =l: ) sm9 i n powers of At, and keeping only terms that survive after take the second 
Euclidean time derivative and coincidence limit, we have 



j&Hng ^ aWZr + ^ A± {3A ± + a) f A _^ , 



At 



6r 2 



-(Ar) 2 + 0(Ar 4 ) , 



(61) 



where we use the short notation = — - — — — . 



The first term on the right hand side of (|6ip provides 



4vrAr (2) ' 

which coincides with the Euclidean bispinor for a flat three-dimensional space-time. 



(62) 
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The VEV of the zero-zero component of the energy-momentum tensor given by (|49p is 
divergent. In order to obtain a finite and well defined result, we should extract its divergent 
part. This can be done in a manifest form by subtracting form the complete Green function the 
usual Euclidean Green function as shown below: 



Ren. 



lim d T TV 

X 1 — >x 



V E {x,x) - V^\x,x) = lim d 2 . Tr V Ren {x, x) . 



So in order to obtain the VEV above we need to develop three different calculations: 



The coincidence limit of the Euclidean time derivative of I{ 



(±)ren , 



I 



(±)sing 
1 



(63) 



T_V2r. 



lim d 2 T l[ ±)ren (z) 



(3^± + 



a 



3r 2 



(64) 



The same calculation as above but now for l[ m (z). In the Appendix |B] we explicitly 
show that 



lim dll[ ±)fin (z) 



1 



being 



/ (±) (y) 





a 

y 



\ f {±) (y) 

dy ^ — 



(65) 



(66) 



The last calculation is for (z). It is easy to show that 



lim dllt\z) 



dy £(±)(y) 



i y 3 y/2 + y* 



(67) 



Finally the renormalized VEV of the zero-zero component of the energy momentum tensor 
can be written in a simpler form, after substituting the expression for A^, using the definition 
for the function f^'(y) above and expressing 2a/3 = 2a dy/y 4 : 



(T$(x))Ren. 



4n 2 aV2r 3 



1 - a 2 



12 7 2 



6 a 



dy 



+ / -|F 1 (y)+ / ^F 2 {y) 



o r 



dy 



where 



Fi{y) 



F 2 (y) 



S(+)(y) + s(-\y) 
SM(y) + S(-\y) 2a 

V^ + y 2 y 



2a 1 
— + — 

y 



i y° 



a 2 - 12<y 2 



(68) 



6 a 



(69) 



The behaviors of the integrals above as functions of a and 7 are exhibited in figure HJ 

Having found Fq the other components of the renormalized VEV of the energy-momentum 
tensor can be expressed in terms of it. Assuming 



<7j 



B I Ren. 



F A 
B 



(70) 



by the diagonal form of the metric tensor, we expect that F$ be a diagonal matrix; moreover, 
the conservation conditions, V ' a{T a ) R en . = 0, provides F$ = —2F£ and (T^)# en . = 0, F® = F r r . 
So we conclude 



F§ = F °diag(l, 1, -2) . 



(71) 
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Figure 1: These graphs provide the behaviors of both integrals in (|68p as functions of a and 
7, admitting these parameters varying in the interval [0.1, 0.9]. The left panel corresponds the 
first integral and the right panel corresponds to the second one. 

3.1.2 N=2 

This case corresponds to four-dimensional one. For this manifold the renormalized VEV of the 
energy-momentum tensor associated with fermionic fields, has been analysed by Linet, cited in 
|25j . and by us in [33]. So we shall not repeat this calculation; however for completeness we 
write down the zero-zero component of this tensor: 

< T °° (X) >- = 28807TW ^ " 1)(17 ° 2 + 7) + 12072( ° 2 " ^ + 1} ] • (72) 

Also, having obtained F® the other components of the renormalized VEV of the energy-momentum 
tensor can be given in terms of it. Assuming a diagonal form for this tensor, a cylindri- 
cal symmetry and invariance under a boost along the x direction, the conservation condition 
V A (T r A ) R en. = 0, provides F£ = -3F r r , and (T^) Ren . = 0, F ° = F r r = F£. So we conclude 

Fb = ^o°diag(l ; 1, 1, -3) . (73) 

3.1.3 N=3 

In braneworld scenario this manifold corresponds to a conical 3— brane having an one-dimensional 
extra space transverse to itH The analysis of the VEV of the energy-momentum tensor asso- 
ciated with scalar field in a higher-dimensional cosmic string space-time and in the presence 
of a magnetic flux has been developed by Guimaraes and Linet in [29]. There, integral rep- 
resentations for massive scalar Green functions have been obtained. Specifically for the four 
dimensional space, the authors calculated the renormalized vacuum expectation value of the 
energy-momentum tensor, (T0 Ren., and show that considering the mass of the field equal to 
zero more suitable expressions are obtained. Moreover, Linet cited in [25], using previous result 
calculated the VEV of the energy-momentum tensor associated with fermionic fields, relating 

6 The five-dimensional cosmic string space-time has been considered in the braneworld theory with metastable 
graviton on the brane in [46] . 
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the parameter 7, the fractional part of the ratio &/&o, with the conicity parameter a. Al- 
though a general formalism has been developed, he only provides the vacuum energy density for 
a massless four-dimensional fermionic field. 

The Dirac matrices for a five-dimensional cosmic string space-time, reads 



r u = 7 (°> r r = f • 7 , r v = — w • 7 , r x = 7 ( 3 W r y = i l5 . (74) 

ar 

Moreover, the Green function associated with this space is the 4x4 diagonal matrix below 



1 



T>e(x,x) 











being 



with 



V {2 \x',x) 



\/2(2tt) 3 

)3/2 dz V 









/<-)(*) 



2=cosh u 



dt 



V cosh t 



dy 



2-1 



(75) 



(76) 



(77) 



' arccoshz 

In the above expression we have introduced the variable t := 2 arcsinh(y/\/2) in the second 
term of the right hand side. As in the analysis developed in subsection 13. 1.1 1 we shall divide the 
integral interval form [yfz — 1, 1] to [1, oo], and define respective integrals by l[ (z) and Jj \z). 
Subtracting and adding into the integrand of l[ (z) the three first terms of the expansion 

S^(y) 



(78) 



where A^ 1 
have: 



^g!Hg(^ and B(±) 



(a 2 -l)(lla 2 +l)+30(8 ± ) 2 ((S ± ) 2 -l) 2 +ma 2 5 ± (l-8 ± ) 



4±\z) 



dy 



+ 2 



/ i 3 T vV - {z 

Vv 2 



S^jy) 
Vi + v 2 ' 



a 

y 



180a 3 



A (±) y _ B (±) V 3 



we 



2-1 

1 ' 



(z-l) 



(±)sing 



The singular contribution to the VEV is contained only in /j =t: ^ m9 (z), which has the form 



(79) 



l[ ±)sm9 (z) 



2a/3 



+ 2V2 



(22-1) 



(80) 



being /3 = arctg ( y f=§ 

Taking the derivative of l| =l: ^ m9 (z) with respect to z and substituting z = coshu = 1 + ^^-0 
we expand the result in powers of Ax which survive after taking the second derivative and the 
coincidence limit. Doing this procedure we get 



dl 



(±)sing 
1 



(*), 



airr 



■'V2 



dz 



2r z 



+ (b&- 
(Ai) 3 V 

(3a + + 5A& 

20r2 



A (±) 



■(Ax) 2 + O(Ax) 4 . 



(81) 



The other two contributions for the Green function are given by: 



7 In this development we have taken first r' = r and tp' — (p. 
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The derivative of l[ in (z) which reads 



dl 



(±)/m. 



where 



dz 



dz 



dy 



- Vy 2 - (z - 1) 



f {±) (y) 



a 

y 



(82) 



(83) 



V2 + ¥ 

As mentioned in Appendix [Bj to take the derivative above it is easier if we introduce a new 
variable b = \fz — 1 into Vy^ in (z) and change y — > by. Doing this procedure we obtain 



dl[ ±)fm (z) 
dz 



s/y 2 



dy d ( f^\y) \ /(*>(!) 
y ) 



b 2 dy 



(84) 



The third contribution comes from the derivative of li; , which easily provides 



d4 ±] (z) 
dz 



dy 



1 (y2 _ fo 2)3/2 



(85) 



The integrals ([84^) and ([85|) have to be evaluated at the point b = 

Now taking into account all these results the Green function can be expressed by 



T> E (x',x) 



where 



\Z2a(27r) 



I KM(Ax) 

K(~\Ax) 

K( + \Ax) 











(86) 



K(-\Ax) j 



dl[ ±)sin9 (z) 



+ 



dl 



(±)fin 
1 



(z) , dl^\z) 



+ 



(87) 

dz dz dz 

The formal expression for the VEV of the zero-zero component of the energy-momentum 
tensor is given by (|49h . However this expression is divergent, and the term responsible for that 
is precisely the first term on the right hand side of (|8ip . which when substituting into into (I75p 
provides 

v ^' x) = ^h?^- (88) 

So, in order to obtain a finite and well defined result to this VEV we should extract its divergent 
part. Here, this also can be done in a manifest form by subtracting form the complete Green 
function the usual Euclidean Green given above: 



(T$(x)) Ren . = lim O 2 Tr V E (x, x 1 ) - vf{x, x') = hm d 2 Tr V Ren .(x,x' 



(89) 



Finally we are in position to write down the the zero-zero component of the renormalized 
VEV of the energy-momentum tensor in a simpler form. Substituting the expressions for 
and BW, using definition for the function f^\y) and some intermediates steps we get: 



(T$(x)) 



Ren. 



[a 



l)(23a 2 - 7) - 120 7 2 (1 - 3a 2 - 2 7 2 ) 



8TT 3 aV2r 5 

Jo y b Ji y b 



240a 2 



(90) 
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where 



F 3 {y) 



+ 



12 7 2 



S(+)(y) + S(~Hy) 2a (1 - a 2 

[(a 2 - l)(17a 2 + 7) + 120 7 2 (1 + a 2 - 2 7 2 )] 
720a 3 ' 
5(+>(2/) + 5(->(y) 2a 

+ y 



(91) 



Here also it is not possible to provide analytical expressions for the integrals above. The numer- 
ical behaviors for both integrals are plotted in figure [2] as functions of the parameters a and 7. 




Figure 2: These graphs exhibit the behaviors of the integrals in (|90p as functions of a and 7, 
admitting these parameters in the interval [0.1, 0.9]. The left panel corresponds the first integral 
and the right panel corresponds to the second one. 

The other components of the energy-momentum tensor can be obtained from the zero-zero 
one, by using the conservation condition (|45p . the absence of trace anomaly (|46p and the boosts 
symmetry in the directions parallel to the string. Writing 

(T£(x)) Ren . = ^ , (92) 

we can easily show that 

Fg = F °diag(l, 1, -4, 1, 1) . (93) 

3.1.4 N=4 

As we have mentioned beforee [3 [9], the gravitational effects associated with strings have been 
considered as responsible for compactification from six to four dimensional space-time. In this 
scenario the bulk is represented by a conical 3— brane transverse to a two dimensional space. 

The manifold analyzed by us here presents a simpler geometrical structure. The bulk corre- 
sponds to an idealized conical 3— brane transverse to a flat two dimensional space. 
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The explicit analysis of the VEV of the energy-momentum tensor developed here is a new 
one. Fortunately, it becomes much simpler when compared with the previous analysis because 
the respective Green function can be expressed in terms of elementary functions. 

In this six-dimensional cosmic string space-time, the Dirac matrices are in the form ([3]), and 
the Green function is expressed in terms of a diagonal 8x8 matrix: 



T> E (x',x) 



( V^\x\x) 













V 






E ( X ' X ) 





(2) 



\ 







Vf{x',x) J 



(94) 



being 



vf{x',x) 



e iN(<p-<p r ) 

5 o s — x 

167r 3 a(rr') 2 sinh u 

cosh uS^ (u) — smhuS'( + \u) 




cosh uS( ' (u) — sinh S'( > ( 



u 



(95) 



with S^(u) having the same functional form as S^(t) given in (|29p . 

As we have already discussed, the calculation of the renormalized VEV of the zero-zero 
component of the energy-momentum tensor can be calculated by applying the second order 
Euclidean time derivative on the renormalized Green function: 



Ren. 



lim cL TV 



X, X 



lim d 2 T Tr V Ren {x,x) , 



(96) 



where T>^\x,x') is promptly obtained from (|94p by taking a = 1 and 7 = 0. After some 
intermediate steps we obtain for the energy density the result below: 



Ctf(*)>T 



1 



120960vr 3 a 6 r 6 ^ 367( ' 1 
+144(357 4 - 497 2 -4)(l 



a 



2\3 



2\2 



a 



+ 12(1897^ + 76) (1 - or) 
) + 1344 7 2 (7 2 -l)( 7 2 -4)] 



(97) 



The other components of the energy-momentum tensor can be computed by using the con- 
servation condition, the absence of trace anomaly and the boost symmetry along the directions 
parallel to the string. So we have: 



(98) 



where 



F£ = F °diag(l, 1, 



-5, 1, 1, 1) 



(99) 



being Fq given in (|9 



3.2 Massive case 

Here we shall analyze the VEV associated with massive field under the circumstance specified in 
12.21 In order to that we shall adopt the bispinor (|42p to calculate the fermionic Green function. 
However, before to start the calculation let us analyse this bispinor in the coincidence limit. 
Taking x' — > x we verify that the result is divergent and that the divergence comes exclusively 
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form the k = component. So, in order to obtain a finite and well denned result we should 
apply some renormalization prescription. This procedure can be applied in a manifest form by 
subtracting from (|42|) its k = component. So the renormalized Green function is given by 



e iNAip I- 1 



(2^)^/2- 



T MN/2 p j^m K m( M P^( I e i(i-i/Jv e -2i. k/q ) • (ioo) 



Before to start the calculation of the renormalized VEV of the energy-momentum tensor, we 
shall analyse the behavior of the renormalized bispinor in the coincidence limit. 
Considering first the upper diagonal component of (jlOOp . we have 

(+) _ M W ^ K N/2 (2MrMk*/q)) 

"Ren\*^) (27r)J V/2+l (2r) JV/2 (sm^TT / q)) N / 2 ' 1 j 

We can see that the T>^ n (x,x) is positive everywhere. Now we would like to analyse this 
expression in two limiting cases: 

• In the limit Mr >> 1/ sm(kn/q), the main contribution comes from k = 1 and k = q — 1 
and the leading term is 



^Ren.{ x i x ) ~ jv+i , , at+1 • (102) 

(47r) 2 (rsin(7r/g)) 2 

• For the massless limit, 

P W (x x) - r (^/ 2 ) V 1 (103) 

For even iV, the summation on the right hand side of (|1U3|) can be developed by using the 
formulas |35| 



M^ir 1 e -2Mrsin(7r/g) 



, , I'/r(x) +N 2 I N (x) t t / \ Q 2 1 / t s 

iV(iv+l) sm (qx) sin (x) 



for the sum 



In particular for a four dimensional space, N = 2, /2(C)) = ^ 3 ^ . So we have 

V$LM = ^- (I 06 ) 

For a six-dimensional space, 1.4(3;) is a long expression obtained by the recurrence relation above, 
however we find /4(C)) = ^ q — ^jf , consequently 

p w fx x) _ (g 2 -i)^ 2 + n) fl07) 

The above results, Eq.s (|106p and (|107p . are analytical functions of q, and by analytical contin- 
uation are valid for all arbitrary values of q. 
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A similar analysis can also developed for the lower diagonal component of (jlUUp . So let 
consider 

V H (x x) _ M»'* £i K N/2 (2Mr S in(kn/ q) ) 

RenX ' ' ( 2vr )^/2+i( 2r )^/2 A. ( s \n{kTi/q)) N / 2 " { ' 

In the two limiting cases above we have 

• For Mr » l/sin(kir/q) we find 

«(-) (r r) „ cos^/gK 2 ^^) 
(47r) 2 (rsm(7r/g)J 2 

• For the massless limit, 

75H _ r(JV/2) e -2ifc7T/q 

^W 3 ^) - (47r) 7V/2+l r 7V sin (fc7r/g))^ " ( ' 

In order to develop the sum in (jllOp for iV even, the relevant formula is 

^ = E-4^=E^§ = M0)-2/ N - 2 (0). (Ill) 
For a four-dimensional space A2 = ^2(0) — 2I (0) = ^ 9 g" 1 ^ — 2(q — 1) = fa -1 ^ 9-5 ) ; consequently 

= '"'jit' 5 ' • ("2) 

For six-dimensional space, If 4 = Zi(0) — 2/2(0) = - ~ 19 ^ and 

P H (x x) - fa 2 - 1 )^- 19 ) (113 ) 

U RenX^)~ 28807r 3 r 4 " 

Now let us proceed the calculation of (T^f (x)) Ren .. Although we shall concentrate the 
analysis of this quantity for a six-dimensional cosmic string spacetime, below we present its 
zero-zero component for a general value of N . By using previous result, 

(T$(x)) Ren . = lim d 2 T Tr V Ren .(x,x') , (114) 

x —*x 

and substituting the Euclidean bispinor expressed by its diagonal form given in terms of (jlOOp . 
after some intermediate steps we obtain 

WW)- = g wtr /, )) »w ' f »^' l °W«)) ■ < 115 > 

From (|115p we can observe that the energy density is non-positive quantity everywhere, vanishing 
only for q = 2. For this case the contributions from the upper and lower component of the 
bispinor (|100p cancel each other. Moreover in the limit Mr >> 1/ sm(k7r/q) the leading order 
is 

(T$(x)) Ren « --^L_. _ ^ x (sin ^ S / ^i 3 )/2 e " 2MrSm(V9) > ( 116 ) 
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with exponentially suppressed behavior. 

In the massless limit case, the energy density became expressed in terms of elementary 
functions, 



r(JV/2 + l) ^4 cos 2 (kn/q) 
' ir(4ir) N / 2 r N + 2 (sm(kir/q)) N + 2 ' 



(117) 



For even N, the summation on the right side of the above equation can be expressed in terms 
of analytical functions of q as shown below: 



T(N/2 + 1) 
'^(4^/2^+2 



(W(0)-/tv(0)) 



(118) 



For N = 4, we easily find from the recurrence formula (|104p . -^(O) = ^(q 2 — l)(2q 4 +23q 2 + 191), 
consequently 



<T °(*)> 



Ren 



(q 2 -l)(q* + q 2 -20) 



(119) 



3780vr 3 r 6 

The above quantity is positive for 1 < q < 2. 

Now going back to the general expression (|115|) . in figure 3 are plotted the behaviors of 
(Tq ) R en ./M N+2 as function of Mr for different values of N and the parameter q. As we have 
discussed, for large value of Mr the energy density tends to zero exponentially. The right panel 
explicitly shows that the modulus of the energy density increases when q becomes larger. 
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Figure 3: These graphs represent the dependences of (Tq) Ren jM N+2 with Mr for different 
values of the parameters N and q. The left panel corresponds to q = 4 and N = 2, 3, 4, and 
the right panel corresponds to N = 4 and q = 3, 4, 5. 

Now after this discussion about the behavior of the energy density for a general value of 
N, let us specialize in a six-dimensional cosmic string spacetime. Below we present the most 
important analytical results: 

Substituting a = - and 7 = Szl into ((97| we re obtain (fTT9)) . 
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Taking advantage of previous result we can promptly write down the zero-zero component 
of the VEV of the energy- momentum tensor. It is 



<2o(s)W 



M 3 l-k cos 2 (kn/q) „ tn% . „ 

z2 3 n rrr K 3 t 2Mr sin (W?)] 



3 r 3 ^ sin 3 {kn/q) 



(120) 



The radial pressure, (TJT), can be calculated as shown below: 

(T rr (x)) Ren . = - lim Tr [T r (d r - d r >)] S F (x, x) 

1 X — *x 

By using fjl2[) to calculate the Feynman propagator, we can write 



(T rr {xy) Ren. 



— lim Tr 

2 x'—>x 



1-cO 



2ar 



(d r — d, 

D E (x, x') 



ar ( °' 



Using eAtp = N + 7, 7 = — a = 1/q and the identities 



(121) 



(122) 



r r r r 



-/ , r r r b 



— s s and r r r y n 6 



2ar 



:i-a)/ 



after a long calculation we find 

(Tr(x)} Ren . -- 



M 3 g — 2 



cos 2 (/c7r/g) 
87r 3 r 3 ^ sin 3 (/c7r/g) 



ir 3 [2Mrsin(fc7r/g)] . 



(123) 



(124) 



From the conservation condition V a(T t ) R en . = 0, the azimuthal pressure, (T<jf), can be 
obtained in terms of the radial one 



CC<»>/?en. = <9 r (r(Tr(x)) Ren .) . 



(125) 



Substituting (I124p into the above expression we get 



COT) 



Ren. 



M 3 q ^cos 2 (kir/q) 



{5i^ 3 [2Mr sin(fc7r/g)] 



87r 3 r 3 sin 3 (A:7r/g) 
+ 2if 2 [2Mrsin(fc7r/g)]Mrsin(fc7r/g')} . 



(126) 



For the pressures along the directions parallel to the string we have (no summation over 





(Ti)Ren. = {T^Ren. , i = 3, 4, 5 



(127) 



Now on basis on the results obtained, we can verify the correct trace of the energy- 
momentum tensor: 



{T£{x)) Ren . 



M 4 " 1 ■■ ■ 



E ^Jnl^l Ki i 2Mr Mkir/q)] = M 2 Tr/;,,, „. ( , . , ) . ( i 2 N ) 



4-7r 3 r 2 s\v?{k-K/q) 
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4 Conclusion and Discussions 



In this paper we have investigated the local one-loop quantum gravity effects associated with 
charged fermionic fields in a higher-dimensional cosmic string spacetime. As the first steps in the 
evaluation of the renormalized VEV of the energy-momentum tensor, in Section [2] we calculated 
the general Euclidean Green function, which is expressed in terms of not-solvable integral for 
a general situation. However, two limiting cases provide more workable expressions which were 
explicitly analyzed: i) the massless case, and u) the massive one with a = |, being q an integer 
number, and 7 = For the latter, we were able to express the Green function in terms of a 
finite sum of modification Bessel functions, K^. 

In Section [3l we considered the VEV of the energy- momentum tensor. For the massless case, 
we calculated, in explicit form, the energy-density for different dimensions. For even dimensions, 
(Tq )nen. is expressed in terms of elementary functions; however for odd dimensions, it is given 
in terms of integrals expressions. For these situations we provided numerical analysis, as shown 
in figures 1 and 2. Being q an integer number and 7 = the corresponding renormalized 
VEV of the energy-momentum tensor is expressed in a closed form for massive fields. To our 
knowledge no such closed expression has been given before in the literature. For an arbitrary 
value for N, we explicitly shown the behaviors of the energy-density for several limiting cases. 
For points near the string, it behaves as l/r N+2 , and for points far away from it, presents an 
exponentially suppressed behavior. Finally, for a six-dimensional space, we calculated completely 
all components of the energy-momentum tensor, and shown they provide the trace identity. 

An interesting point which deserves to be mentioned is that the effects on the renormalized 
VEV of the energy-momentum tensor due to the conical structure of the spacetime and the 
magnetic interaction, may cancel each other. This was explicitly observed in Subsection 13.21 for 
q = 2 and 7 = 7, and can be also verified in the results found for the energy-densities given 
in ([72|) and ([97|) . In fact, for these specials values, and taking Ay? = 0, the upper and lower 
components of (f3T|) can be expressed as 



V E \x',x) 



1 cosh u 

47r 2 r 2 sinh 2 u 
1 1 

47r 2 r 2 sinh 2 u 



(129) 



with cosh(it) given by ([25]) . which coincides with the respective quantities found by image 
method, 



-D { +\x\x) 



V E \x',x) 



1 


" 1 


1 " 

+ — 
Pl. 


4^2 


A 


1 


' 1 


1 " 


4^2 


A 


~7i. 
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being po and p\ given in (|43p . In this situation the contributions on the renormalized VEV of the 
energy-momentum tensor due to the upper and lower components of the Feynman propagator 
cancel each other. 

The renormalization procedure adopted in this paper to calculate the vacuum expectation 
value of the fermionic energy-momentum tensor was the point-splitting one. In this sense we had 
analysed the behavior of this quantity at the coincidence limit and extracted all the divergences 
in manifest form by using the Hadamard function]! It has been proved that this procedure 
provides finite and well defined results [IT]- An alternative regularization procedure can also 
be used. The local (^—function procedure leads essentially to the same conclusions as given 



Because the cosmic string spacetime is locally flat the Hadamard function coincides with the usual Green one. 
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by the point-splitting one [IS]. The DeWitt-Schwinger expansion of the Green function for an 
n— dimensional curved spacetime is given in |41|I44|. So, by using the dimensional regularization, 
it follows from the general expression that for the cosmic string spacetime this expansion contains 
a single term only, 

yDStj „\ _ in ( z V~ n/2 o-(2) 



G ' (4^ ' 131) 

where is the Hankel function of the second kind, z 2 = —2m 2 a being a the one-half of the 
square of the geodesic distance. The above function coincides with the Minkowski Green func- 
tion. The subtracted contribution due to this function renormalizes the cosmological constant. 
All these regularization procedures lead to the same finite result for renormalized quantities. 
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A The generalization formula of the generating function of mod- 
ified Bessel Function I v 

The method of images applied in (|2.2.ip to calculate the Green function associated with massive 
fermionic fields in a conical space-time and in the presence of a magnetic flux running along the 
string, is an application of the generalization formula of the generating function for the modified 
Bessel functions of integer order. It well known that 



oo 

,x/2(t+l/t) 



Y w* n • ( 132 ) 



Taking t = e lip we have 



'OO 

x COS tf> 



Y u^y^ . (133) 



Making the change 99 — > ip + in (|133p . and summing over k from until q—1, we may write 

q—1 00 q—1 

fc=0 n=-oo k=0 



Now using the identity 



we obtain 



q—1 00 

Y^ nk/q = q Y (135) 

k=0 m=-oo 



q—1 00 

Y^e xcos ^ +27rk/q) =q Y I mq (x)e im ^ . (136) 

k=0 m=-oo 
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B Explicit calculation of F ° 

In order to find the expression (|65p for the limit case analysed, our first steps is to define a new 
variable 6 = \fz — 1 = In terms of this variable we may express l[^^ m (z) as: 

l[ ±)fm (z)=2 (137) 
Jb \Jy l - b 1 

being f^\y) given by (|66p . Using the fact that d 2 = (l/2r 2 )<9 2 we have 

*f^w = ^f Jh^- (138) 

Because the integrand is divergent at the point y = b, we have to change the variable y — > by 
before applying the Leibnitz formula for derivative of an integral. So the second derivative has 
the form below, 

dPJ b y^V 1 {V) PJb yW^ J [V> b 2 VT=V 62(1-62)3/2 ' ^ 
which can be written as 

*2 7 fe - &2 y 6 ^2 _ 52 1 y y 2 J 



f( ± )(1) /( ±) (1)6 : 



2 



V 7 !^ 2 " (1 - 62)3/2 

Taking the limit 6 — > we obtain 



(140) 



*r * f*^uq*. (14 i) 

«62 J 6 ^7/2 _ 52 J 



We can see that the integral above is finite because for small value of y, (y) 
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